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ABSTRACT 

We combine HST/NICMOS imaging photometry of the HR 4796A disk at previously unobserved wavelengths 
between 1.71-2.22/im with reprocessed archival observations to produce a measure of the dust's scattering 
efficiency as a function of wavelength. The spectrum of the dust, synthesized from the seven photometric 
measures, is characterized by a steep red slope increasing from 0.5 /im to 1.6 /mi followed by a flattening of 
the spectrum at wavelengths > 1 .6 /im. We fit the spectrum with a model population of dust grains made of 
tholins, materials comprised of complex organic materials seen throughout the outer parts of our Solar System. 
The presence of organic material around a star that may be in the later stages of giant planet formation implies 
that the basic building blocks for life may be common in planetary systems. 

Subject headings: stars:individual(HR 4796A)-circumstellar disks 



1. INTRODUCTION 

HR 4796A is a young (8 Myr) star with a circumstel- 
lar debris disk that has been imaged at several wavelengths 
(Jayawardhana et al. 1998; Koerner et al. 1998; Schneider 
et al. 1999, 2008). Dust in the disk scatters starlight at visible 
to near-IR wavelengths and emits thermally in the infrared. 
Both scattered and thermal emission are spatially coincident 
originating in a < 17 AU wide annulus (Jayawardhana et al. 
1998; Schneider et al. 1999; Koerner et al. 1998; Schneider 
et al. 2008) 70 AU from the central star. Estimates of the 
mass of material present in the disk range from 10 26 - 10 27 g 
(Jura et al. 1993, 1995). The constituent grains assumed to be 
responsible for the spatially resolved scattered light and ther- 
mal emission have radii on the order of a few /im although 
unresolved emission at infrared and sub-mm wavelengths re- 
quire larger grains a few tens of /im in size to be present 
(Augereau et al. 1999; Li & Lunine 2003; Wahhaj et al. 2005; 
Kessler-Silacci et al. 2005). The smaller grains will escape 
from the HR 4796A system on short timescales implying a 
resevoir of colliding bodies that exists to fuel the dust pro- 
duction. The characteristics of the HR 4796A disk make it a 
younger, dustier analogue to the Kuiper Belt within the So- 
lar System. Studying HR 4796A provides a probe of planet 
formation processes around a star with a mass and effective 
temperature significantly different than our own. 

The ring structure of the disk was inferred by thermal emis- 
sion well fit by a single blackbody temperature (Jura et al. 
1993, 1998). It was recognized that the disk could contain 
dust that was not simply made of silicates, but available ob- 
servations could not usefully constrain the true composition 
of the dust grains (Augereau et al. 1999; Wahhaj et al. 2005). 
Mid-IR Keck LWS spectra from 8-13/im similarly showed a 
lack of observable silicate lines or other identifying features 
(Kessler-Silacci et al. 2005). Near-IR scattered light images 
of the circumstellar dust (Schneider et al. 1999) suggested a 
red spectral slope at these wavelengths, hinting at a compo- 
sition similar to that seen around Solar System bodies in the 
Kuiper Belt. 

2. OBSERVATIONS 
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In order to measure the reflectivity of the dust, we observed 
HR 4796A with the Hubble Space Telescope's (HST) Near- 
Infrared and Multi-Object Spectrometer (NICMOS) (Thomp- 
son et al. 1998) using the coronagraph in Camera 2. We 
used four medium bandwidth filters (AA/A ~ 10%), F171M, 
F180M, F204M, and F222M, with A c = 1.71, 1.80, 2.04, and 
2.22 /im. 

All new observations were taken on 09 May 2005 for 
HR 4796A and HR 4748, a PSF reference. The observations 
include direct images of both stars with short exposures for 
stellar photometry, as well as longer exposures for corona- 
graphic high contrast imaging. The instrumentally calibrated 
and reduced images were created from the raw NICMOS 
multiaccum exposures following the processing methodolgy 
described by §3 of Schneider et al. (2005) and references 
therein. Additionally, we took archival raw images in the 
STIS 50CCD (A=0.585), NICMOS F110W (A=l.l/Ltm), and 
F160W filters (A=1.6/im) and re-reduced the images. 

For each of the medium bandwidth filters we directly mea- 
sured HR 4796A and HR 4748's flux densities with their 
uncertainties to obtain wavelength dependent brightness ra- 
tios for PSF subtraction. Due to the (detector saturating) 
brightness, unocculted stellar images were not available in the 
broader passbands. Thus, we estimated the 50CCD, F110W, 
and F160W photospheric flux densities with synthetic pho- 
tometry using the CALCPHOT task of the SYNPHOT 1 pack- 
age in conjunction with Kurucz models tabulated in the SYN- 
PHOT library. We used a solar metallicity, logg=4.5 model 
of a 9250 K star, using a luminosity of 21.1 L Q from a best 
fit match to HR 4796A based on its parallax from the Hip- 
parcos satellite, its Tycho-2 B and V photometry, as well 
as its 2MASS J, H, and K s photometry. For HR 4748 
we used data from the same sources to derive a 121 L Q , 
logg=4.5 model with T<,/y=10250 K. We compared the pre- 
dicted F171M, F180M, F204M, and F222M photometry with 
the SYNPHOT values-they agreed to within a few percent 
of our measured flux densities. We fit the long wavelength 
emission from the disk to a single temperature blackbody and 
found T e// =78± 1 K and L /R /L*=4.8±0.1 x 10" 3 , a useful sur- 
rogate for the disk optical depth. Since <C 1, it is rea- 
sonable to assume that the disk is optically thin. However, if 
the ring had a small enough scale height, the disk's radial op- 

1 http://www.stsci.edu/resources/software_hardware/stsdas 
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FIG. 1 . — False color image of the HR 4796A disk with three of the pass- 
bands imaged. Blue corresponds to the 50CCD image, green is the F110W 
image, and red is the F171M image. A slight color difference between the 
two lobes of the disk is seen. The central region where PSF subtraction errors 
dominate is masked out. 



tical depth would approach unity. We discuss the implications 
of a radially optically thick disk on our proposed model in §3. 

In order to determine the best subtraction we minimized a 
chi-square metric on a region of the diffraction spikes for the 
target star. We assume that good subtraction of the diffraction 
spikes corresponds to the best subtraction of the PSF within 
the region of interest (Schneider et al. 2001). We searched 
within 1-er of our scaling uncertainty and iteratively searched 
to within 0. 1 pixels of the best offsets to get the optimal sub- 
traction between target star and PSF. 

To quantify the systematic effects on the photometry, we 
repeated the subtractions varying the PSF scalings and offsets 
by ±1 (7 from the minimum chi-square solution found above. 
Using an elliptical photometric aperture matched to the size 
and shape of the disk, we found the standard deviation in the 
disk flux densities from this suite of subtractions and used this 
as a systematic uncertainty per pixel that was then propagated 
into the disk photometry. 

We observed HR 4796A in each medium band filter at two 
celestial orientations differing by 26°. This is essentially an 
azimuthal dither that allows true objects that rotate with the 
field of view to be distinguished from rotationally invariant 
instrumental artifacts. The two differentially-oriented PSF- 
subtracted images in each filter band were geometrically rec- 
tified and corrected for the linear optical distortion {X sca i e = 
75.95 mas pixel" 1 , Y jefl / e =75.42 mas pixel" 1 ) at the Camera 
2 coronagraphic focus. The geometrically corrected image 
pairs were rotated about the location of the occulted target 
(as determined through the target acquisition process), then 
co-aligned to higher precision using the image centroids of a 
nearby field star common in both images, and, finally, com- 
bined by averaging. 

To determine the brightness of the disk from the PSF sub- 
tracted images, we chose an elliptical aperture with a major- 
to-minor axis ratio equal to the measured inclination of the 
disk of 14.1° from edge-on (Schneider et al. 2008). We used 
an inner major axis boundary of 0"8, and an outer major axis 
boundary of l."5. At r <0"8, systematic errors dominate the 
image, which restricts useful photometry of the disk to within 
±20° of each ansa. Beyond 1"5, disk flux is below 5-10% 
of the peak. We used an elliptical annulus between 1 "9 and 



2"5 with the same azimuthal coverage to determine the back- 
ground level of the image. Table 1 gives the total measured 
flux densities of the disk for the four filters. 

There were three differences with our re-reduction and sub- 
traction of the PSF star for the Cycle 7 F110W and F160W 
data. We directly calculated pixel centers of the stars behind 
the coronagraphic hole using the acquisition images and slews 
reported in the data file headers. A flight software problem 
caused the relevant header keywords of the acquisition image 
centroid to be incorrectly populated. The geometric distortion 
of Camera 2 is time dependent, so we used X JtY ,/ f =75.87 mas 
pixel" 1 and Y jefl ; e =75.19 mas pixel" 1 , appropriate for these 
observations. 

We measured disk asymmetries from our new and re- 
prossessed images, as has been done previously for 50CCD 
and Fl 10W (Schneider et al. 2008; Augereau et al. 1999). The 
disk appears to be asymmetric about the major axis, in a man- 
ner consistent with forward scattering dust grains. We quan- 
tify this asymmetry in the context of the Henyey-Greenstein 
g parameter. Additionally, the northeastern (NE) side of the 
disk is brighter than the southwestern (SW) side. We present 
both sets of measurements in Table 1. There is a marginal 
trend to a larger g- value that peaks at 1.71/im. The asymme- 
try between the two sides of the disk is particularly strong at 
1.71 and 1.8/im 

Figure 1 shows a three-color image of HR 4796A using in- 
formation from the STIS, Fl 10W, and F171M filters to high- 
light visually the spectrum of the disk and any spatial varia- 
tions in the spectrum. 

To determine the assumed forward scattering asymmetry of 
the ring, we measured the total flux in the four readily visi- 
ble quadrants (between 5 and 15° from the major axis) of the 
disk between 0."8 and 1 "5, and took a ratio between the front 
(brighter) and back (dimmer) parts of the NE and SW ansae 
of the disk. We calculated the same ratios using a Henyey- 
Greenstein phase function (Henyey & Greenstein 1941) for 
the disk inclination to determine a best fit value for the asym- 
metry g at each wavelength. To measure the brightness asym- 
metry between the two lobes, we took the total flux in the NE 
and SW lobes and measured their ratios (±20° from the ma- 
jor axis) between 0''8 and 1"5, using the same background 
annulus. 

It is useful to determine how the observed surface bright- 
ness in an optically thin disk relates to dust properties. The 
average surface brightness (F v ^isk/Adisk) of an optically thin 
disk is the product of the total scattering cross section of the 
dust integrated over the volume of the disk multiplied by the 
flux density of the star (F^ + ) and the phase function of the 
dust $(#,#) (Augereau & Beust 2006; Schneider et al. 2006). 
Since our observations were close to the ansae, we assumed 
that the phase function of the dust could be approximated 
by the Henyey-Greenstein phase function at a scattering an- 
gle of ~ 90° and assuming the median g value for all of the 
wavelengths. The measured asymmetries may be used to cal- 
culate a geometrical factor relating the total disk flux to the 
observed flux close to the ansae. The geometrical factor was 
calculated by using an azimuthally averaged radial cut of both 
lobes of the disk (±20° from the major axis) and by generat- 
ing a model disk from that measurement. Additionally, we 
multiplied the model disk by a sine function to approximate 
the observed NE to SW asymmetry. We measured the ratio of 
the total flux in the model disk to that of the apertures we used 
to determine the geometrical factor at each wavelength. 
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Table 1 . Properties of the HR 4796A disk 



Filter 


Scaling 2 


Flux Density 


Total Disk Flux Density 


NE:SW ratio 


g 








yuuy ) 






STIS 50CCD 


0.705 


5.5±0.2 


9.4±0.8 


0.74±0.07 


0.16±0.06 


F110W 


0.78±0.03 


5.2±0.3 


8.1±0.6 


0.9±0.08 


0.06±0.07 


F160W 


_b 


3.7±0.5 


5.7±0.9 


1.0±0.1 


0.0±0.06 


F171M 


0.819±0.007 


2.8±0.1 


4.8±0.5 


0.62±0.08 


0.3±0.1 


F180M 


0.79±0.01 


2.9±0.1 


4.7±0.3 


0.71±0.07 


0.2±0.1 


F204M 


0.79±0.02 


2.5±0.1 


4.2±0.2 


1.02±0.09 


0.2±0.1 


F222M 


0.79±0.03 


2.1±0.1 


3.5±0.2 


1.0±0.1 


0.1 ±0.1 



"Ratio of HR 4796A to HR 4748 

b No contemporaneous PSF reference, used five different references with different scalings 



3. HR4796A'S DISK SPECTRUM 

A combination of the surface brightness and Ljr/L* in con- 
junction with Mie scattering models of a given dust distribu- 
tion and composition can self-consistently provide a measure 
of the scattering efficiency {(Q sca )) of the HR 4796 A disk: 



'.disk 



2\1.5 



<a 2 >A disk F v ^ (l-g 2 ) 



J {<Tqbs)F v ^dV 

J F V{k dv 



Lir 
(1) 

where A c n s k is the area of the disk, g is the asymmetry pa- 
rameter for the dust, tt (a 2 ) is the effective grain cross section 
given a power law dust distribution with a minimum and max- 
imum grain size, and (a_&s) is the effective absorption cross 
section per grain of the dust. 

Figure 2 shows F v ^iskjF v ^ spanning from the visible to the 
near-IR. The red slope of the disk is reminiscent of organic 
materials on the surface of Centaur 5145 Pholus but not of 
other expected components of planetesimals or common red 
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FIG. 2.— F V:ibk /F Vl + for HR 4796A. Horizontal bars indicate the FWHM 
of each bandpass. The red curve and symbols represent 5145 Pholus. The 
other curves represent grain models for candidate materials with an n~ 3 ' 5 size 
distribution with c/„„„=3^tm and a„ MV =1000/xm. Symbols represent the mod- 
els integrated over each bandpass. Each curve is normalized to the Fl 10W 
data of HR 4796A and offset for clarity. 




materials such as water ice, astronomical silicates, hematite 2 
or olivine reddened by UV ablation (Cruikshank et al. 1998 
Poulet et al. 2002; Laor & Draine 1993; Pollack et al. 1994 
Brunetto et al. 2007), also plotted in Figure 2. The scatter- 
ing properties of organic materials formed in the outer So- 
lar System are not well known, but optical constants exist for 
tholins, complex organic polymers that reproduce the spectral 
characteristics of organic materials in the atmosphere of Ti- 
tan and which have been synthesized on Earth (Khare et al. 
1984). Tholins have been used to successfully reproduce the 
red spectra of icy Solar System bodies (i.e. Cruikshank et al. 
1998). Similar bodies may be colliding to create the dust 
around HR 4796A. 

The top panel in Figure 3 shows a comparison of the disk 
spectrum with a best fit grain model of Titan tholins only with 
two free parameters: a mi „=1.4/im a maA =980^im and a fixed 
power law index of k=-3.5 (Dohnanyi 1969), with a reduced 
chi-squared value of the fit of 1.2. The chi-squared contours 
imply that minimum grain sizes of between 1 /im and 5/im are 
possible within the 98% confidence contour. The spectrum is 
weakly dependent on the maximum grain size and so is not a 
meaningful constraint on the maximum size of the dust. This 
is true for all the models shown. 

Given that the HR 4796A disk's thermal emission can 
be modeled with a 78 K blackbody, it is conceivable that 
methane, water, or some other ice may be present in the disk 
(but see Grigorieva et al. 2007). The spectral shape we see 
may be due to the presence of tholin inclusions in a larger ice 
or silicate matrix, since pure ice and silicate models do not 
possess the red slope seen in the data. Using effective optical 
constants for water ice with tholin inclusions (Bohren & Huff- 
man 1983) and forsterite with tholin inclusions we compare 
the resulting {Q sca ) with our data in Figure 4. The reduced 
chi-squared value for these two mixtures is 6.8 and 2.0, re- 
spectively. These values make an ice+tholin mixture unlikely, 
but do not rule out the presence of tholin and silicate mixtures. 

We note that the general red spectral shape of the disk is 
very similar to 5145 Pholus, one of the redder objects in the 
Solar System, which is believed to be comprised of water ice 
with tholin inclusions mixed with olivine and methanol ice 
(Cruikshank et al. 1998; Poulet et al. 2002). Additionally, we 
confirm that the minimum grain size in our models for the 
scattered light is comparable in size to the minimum grain 
size predicted by a detailed fit of the long wavelength SED 
(Wahhaj et al. 2005). Tholins provide a useful proxy for the 

Optical constants from http://www.astro.uni- 

jena.de/Laboratory/OCDB/data/oxide/hem_interp.dat 
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FIG. 3. — Comparison of three models of the dust disk around HR 4796A- 
(top) tholins with a,„,„=1.4^tm, a,„ nv =98()/im, (middle) water ice with tholin 
inclusions that comprise 10% of the volume with a„„„=12 /im,a„ MV =l()0^m, 
(bottom) forsterite with tholin inclusions that comprise 10% of the volume 
with a m ,„=12/jm,a„ w ,=l iOfim. 



dust present in HR 4796 A, but may not be the exact substance 
observed, given the large number of different compounds that 
make up tholins. 

If, as discussed in §2 the disk is marginally optically thick in 
the radial direction, the overall choice of composition would 
remain unchanged, though the details of grain size and dis- 
tribution may be different. In a multiple scattering case, the 
surface brightness of the disk would depend on the extinction 
efficiency of the dust and g in addition to (Q sca )- The ex- 
tinction efficiency and g don't appreciably vary as a function 
of wavelength betwee 0.5-2^m for tholin grains. Therefore, 
even in the case of multiple scattering, the scattering spectrum 
would depend solely on (Q sca ). Given that the total mass of 
dust is uncertain by an order of magnitude and the scale height 
of the disk is unknown, the optically thin assumption is justi- 
fied. 

The calculations we performed also make strong predic- 
tions on the forward scattering behavior of the disk grains. 
The scattering phase asymmetry parameters predicted by our 
models for tholin dust are much higher (g ^0.8) than actually 
observed (See Table 1). For spherical grains, g is strongly 
dependent on grain size relative to observed wavelength, and 
more weakly dependent on composition. The measured asym- 
metries are much smaller than predicted by our grain model. 
It is possible that the shape and porosity of the grains make 



Augereau, J.-C. & Beust, H. 2006, A&A, 455, 987 

Augereau, J. C, Lagrange, A. M., Mouillet, D., Papaloizou, J. C. B., & 

Grorod, P. A. 1999, A&A, 348, 557 
Bertini, I., Thomas, N., & Barbieri, C. 2007, A&A, 461, 351 
Bohren, C. F. & Huffman, D. R. 1983, Absorption and scattering of light by 

small particles (New York: Wiley, 1983) 
Brunetto, R., Roush, T. L., Marra, A. C, & Orofino, V. 2007, Icarus, 191, 381 
Cruikshank, D. P., et al. 1998, Icarus, 135, 389 
Dohnanyi, J. W. 1969, J. Geophys. Res., 74, 2531 

Grigorieva, A., Thebault, P., Artymowicz, P., & Brandeker, A. 2007, ArXiv 
e-prints, 709 

Henyey, L. G. & Greenstein, J. L. 1941, ApJ, 93, 70 

Jayawardhana, R., Fisher, S., Hartmann, L., Telesco, C, Pina, R., & Fazio, 
G. 1998, ApJ, 503, L79+ 



the assumption of spherical symmetry for Mie scattering in- 
correct and thus markedly changes the phase function of the 
dust. Porous spherical grains or fractal aggregates tend to 
have higher g values (Wolff et al. 1998; Bertini et al. 2007), 
however, so it is not clear what mechanism is responsible for 
the lower values of g observed. Other disks seem to show low 
asymmetry (Schneider et al. 2006). 

4. CONCLUSIONS 

We have measured a visible to near-IR photometric re- 
flectance spectrum of the dust ring present around HR 4796A. 
The scattered light from the ring is suggestive of a grain popu- 
lation dominated by 1.4^m-radius grains of complex organic 
material similar in optical properties to the material observed 
in the atmosphere of Titan and the surfaces of icy bodies in 
the Solar System. The model grains fit the scattered light data 
but do not fit the observed azimuthal asymmetries in the ob- 
served brightness of the dust ring, suggesting that the grains 
are not spherical. 

The library of optical constants for material is not wholly 
complete, but tholins provide the best fit of the current se- 
lection of materials as illustrated in Figure 2. Many miner- 
als and ices can be summarily rejected because they possess 
relatively neutral scattering efficiencies from 0.5-1.5^m. UV 
ablated olivine has been used to explain particularly red as- 
teroids due to space weathering (Brunetto et al. 2007) and 
hematite is found on the surface of Mars, but these materi- 
als cannot reproduce the observed slope we see in the regime 
of Mie scattering. Thus, the presence of organic material is 
the most plausible explanation for the observations. 

Longer wavelength scattered light observations will further 
constrain the grain models we have used, particularly around 
3.8-4/im where a large absorption feature is seen for differ- 
ent grain sizes of tholins. This would help to directly confirm 
whether Titan tholins are an adequate proxy for the material in 
orbit around HR 4796A. Additionally, measuring the optical 
properties of organic materials in meteorites and from sam- 
ples of the Stardust mission will provide further tests of our 
model grains. 
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a grant from STScI. AJW also acknowledges support from the 
NASA Astrobiology Institute. J.D wishes to thank Jeff Cuzzi 
for insightful conversations on Mie scattering and Alexandra 
Surcel for a careful reading of the manuscript. 
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